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The papain-like protease (PL pro ) of Middle-East respiratory syndrome coronavirus (MERS-CoV) has 
proteolytic, deubiquitinating, and delSGylating activities. The latter two are involved in the 
suppression of the antiviral innate immune response of the host cell. To contribute to an 
understanding of this process, we present here the X-ray crystal structure of a complex between 
MERS-CoV PL pro and human ubiquitin (Ub) that is devoid of any covalent linkage between the two 
proteins. Five regions of the PL pro bind to two areas of the Ub. The C-terminal five residues of Ub, 
RLRGG, are similar to the P5-P1 residues of the polyprotein substrates of the PL pro and are 
responsible for the major part of the interaction between the two macromolecules. Through site- 
directed mutagenesis, we demonstrate that conserved Asp165 and non-conserved Asp164 are 
important for the catalytic activities of MERS-CoV PL pro . The enzyme appears not to be optimized 
for catalytic efficiency; thus, replacement of Phe269 by Tyr leads to increased peptidolytic and 
deubiquitinating activities. Ubiquitin binding by MERS-CoV PL pro involves remarkable differences 
compared to the corresponding complex with SARS-CoV PL pro . The structure and the mutational 
study help understand common and unique features of the deubiquitinating activity of MERS-CoV 
PL pro . 
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INTRODUCTION 

To date, six coronaviruses infecting humans have been 
characterized. Infections with human coronaviruses 
(HCoVs) 229E (Hamre & Procknow, 1966), OC43 
(McIntosh et ah, 1967), NL63 (van der Hoek et ah, 
2004), and HKU1 (Woo et ah, 2005) cause relatively 
mild symptoms in most cases, whereas severe acute res- 
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piratory syndrome coronavirus (SARS-CoV; Drosten et 
ah, 2003; Ksiazek et ah, 2003; Kuiken et ah, 2003; Peiris 
et ah, 2003) and Middle-East respiratory syndrome 
coronavirus (MERS-CoV; Zaki et ah, 2012) are connec¬ 
ted with severe respiratory-tract infection and, in particu¬ 
lar in case of MERS-CoV, acute renal failure (Eckerle et 
ah, 2013), leading to high case-fatality rates of-10 and 
-35%, respectively. In spite of 13 years of research on 
SARS-CoV (Hilgenfeld & Peiris, 2013), no approved 
drugs or vaccines are available for the treatment or pre¬ 
vention of coronavirus infection (Wang et al., 2016). 
This is mainly due to the fact that although these emer¬ 
ging viruses have devastating effects on those infected, 
the absolute numbers of cases (-8000 for SARS, 1733 so 
far for MERS; (http://www.who.int)) imply that the de- 
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velopment of specific antivirals is very likely not com¬ 
mercially viable. On the other hand, the global risk posed 
by MERS-CoV must not be underestimated. Since its 
discovery in September 2012, the number of MERS 
cases reported has been rising steadily, with some inter¬ 
mittent peaks connected to hospital outbreaks in Saudi 
Arabia (Assiri et al., 2013). Man-to-man transmission of 
MERS-CoV has also been impressively demonstrated by 
the recent outbreak of MERS in South Korea, which was 
traced back to a single traveller from the Arab peninsula 
(Butler, 2015). Therefore, it is imperative that academic 
laboratories help increase the preparedness against a pos¬ 
sible MERS-CoV pandemic by characterizing antiviral 
drug targets and by identifying lead compounds interfer¬ 
ing with them. 

In order to successfully infect humans, a virus has to 
meet at least two conditions: 1), it should maintain a suf¬ 
ficiently correct replication of its genetic material; 2), it 
should inhibit the host antiviral response. The papain¬ 
like protease (PL pro ) of MERS-CoV (or SARS-CoV) is 
involved in both of these tasks (Yang et al., 2013; Bar- 
retto et al, 2005). The PL pi ° is a domain located in the 
middle part of the largest non-structural protein, Nsp3, of 
MERS-CoV (or SARS-CoV). It is responsible for releas¬ 
ing Nspl, Nsp2, and Nsp3 from the polyproteins la 
(ppla) and lab (pplab), an essential step of replication 
(Harcourt et al., 2004). Like its SARS-CoV counterpart, 
the MERS-CoV PL pi ° also has deubiquitinating (DUB) 
and delSGylating activities in vivo as well as in vitro 
(Yang et al., 2013; Mielech et al., 2014; Lei et al., 2014; 
Baez-Santos et al., 2014b). K48- and K63-polyubiquitin 
poly (Ub) and ISG15 (interferon-stimulated gene 15) - 
conjugated targets are usually involved in host innate im¬ 
mune regulation (Liu et al., 2005; Maringer & Fernan- 
dez-Sesma, 2014). The PL pi ° has the ability to digest 
K48- and K63-linked polyUb chains and to remove 
ISG15 from ISG15-linked proteins (Baez-Santos et al., 
2014b), thereby interrupting the signalling pathways 
leading to the innate immune response. Thus, the PL pi ° 
can block the activation of IFN regulatory factor 3 
(IRF3) (Yang et al., 2013) and subsequently the produc¬ 
tion of interferon p (IFNP) (Mielech et al., 2014). Inter¬ 
estingly, MERS-CoV PL pi ° shows a similar cleavage rate 
for K48- and K63-linked polyUb chains, while the 
SARS-CoV enzyme prefers K48- over K63-linked 
chains (Baez-Santos et al., 2014b). The former enzyme 
degrades a polyUb chain by removing mono-Ubs, where¬ 
as the latter cleaves di-Ub units off the polyUb chain 
(Bekes et al., 2015). 

We have reported the first crystal structure of the 
MERS-CoV PL pi ° (Lei et al., 2014). Later, two other 
groups also described the structure of this enzyme 
(Bailey-Elkin et al., 2014; Lee et al., 2015). The struc¬ 
ture of PL pi ° can be divided into two parts: a ubiquitin- 


like (Ubl) domain and a catalytic domain with thumb, 
palm, and fingers subdomains. The overall fold of 
MERS-CoV PL pi ° is not only similar to that of SARS- 
CoV PL pi °, but also to that of several human ubiquitin- 
specific proteases (USPs) (Hu et al., 2005). In 2014, the 
X-ray structure of the complex of SARS-CoV PL pi ° with 
ubiquitin has been reported (Chou et al., 2014; Ratia et 
al., 2014). Several key residues (such as Glul68 or 
Tyr265) of SARS-CoV PL pi ° that are important for ubi¬ 
quitin recognition (Chou et al., 2014; Ratia et al., 2014), 
are not conserved in MERS-CoV PL pi °. Bailey-Elkin et 
al. (2014) described the structure of an artificially linked, 
covalent complex between ubiquitin and MERS-CoV 
PL pro . Here, we present the crystal structure of a non-co- 
valent complex between the two proteins and a mutation¬ 
al study of the interactions involved. For these studies, 
we used the CyslllSer active-site variant of MERS- 
CoV PL pro . 

MATERIALS AND METHODS 

Recombinant production of MERS-CoV PL pro 
and its variants 

The PL pro of MERS-CoV (strain 2c EMC/2012 ; Gen- 
Bank no. AFV09327.1) contains 320 residues, from 
Glnl482 to Aspl801 of ppla/lab. For simplification, Glnl482 
was renumbered into Glnl here. The DNA plasmid cod¬ 
ing for MERS-CoV PL pro was produced earlier (Lei et 
al., 2014). 

The MERS-CoV PL pro Cl 1 IS, D164E, D164A, D165E, 
D165A, and F269Y variants were produced using the 
same strategy that we described before (Lei et al., 2014). 
All primers for these variants are listed in Supplemental 
Table S . All DNA plasmids coding for the altered PL pi ° 
were verified by sequencing. 

Genes coding for wild-type (WT) MERS-CoV PL pi ° 
and for its variants were expressed and the correspond¬ 
ing proteins were purified according to our previous de¬ 
scription (Lei et al., 2014). 

Recombinant production of SARS-CoV PL pro 

The PL pr °of SARS-CoV (Strain TOR2 ; GenBank no. 
AY274119.3) comprises 319 amino-acid residues, cor¬ 
responding to Glul541 to Tyrl859 of ppla/lab. A gene 
coding for the SARS-CoV PL pi °was amplified by PCR 
with the following two primers 5-CTAGCTAGCGAGGTTA 
AGACTATAAAAGTGTTC-3' (forward) and 5'- 
CCGCTCGAGTTAATACGACACAGGCTTGATGGTTG 
TAG-3' (reverse). The PCR product was digested by Nhe 
I and Xho I, then was ligated into the pET-28a plasmid 
(Novagen). The recombinant plasmid DNA was verified 
by sequencing. Expression of the gene construct coding 
for SARS-CoV PL pi ° and purification of the protein were 
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performed according to the procedure described for 
MERS-CoV PL pro (Lei et al, 2014). 

Crystallization of MERS-CoV PL pro with ubiquitin 

Purified MERS-CoV PL pi ° (WT) and a variant that had 
the active-site Cysl 11 replaced by Ser (PL pi ° (Cl 1 IS)) 
were both concentrated to -24 mg/mL in 20 mmol/L 
Tris-HCl, 150 mmol/L NaCl, pH 8.8, 10 mmol/L P-mer- 
captoethanol (BME). Human ubiquitin (BostonBiochem) 
was dissolved to 6 mg/mL in 20 mmol/L Tris-HCl, 150 
mmol/L NaCl, pH 8.8. 500 pL PL pro (WT) or PL pro (Cl 1 IS) 
were mixed with 500 pL ubiquitin (-1:1 molar ratio) at 4 
°C overnight. The complex of PL pi ° (WT)-Ub or PL pi ° 
(Cl 1 lS)-Ub was purified by gel filtration (HiLoad™ 
16/60 S200 column, GE Healthcare) in 20 mmol/L Tris- 
HCl, 150 mmol/L NaCl, pH 8.8 the next day. The final 
concentration of PL pro (WT)-Ub or PL pro (Cl 1 lS)-Ub 
was -12 mg/mL. The two complexes were crystallized 
using the sitting-drop method and a Phoenix crystalliza¬ 
tion robot (Art Robbins) at 18 °C. 0.25 pL of protein and 
0.25 pL of reservoir were mixed and equilibrated against 
75 pL reservoir. Screening kits Index™, SaltRx™, PEG 
Rx M 1 & 2, PEG/Ion™ 1 & 2 (Hampton Research), and 
Structure Screen 1 MD1-01, Structure Screen 2 MD1-02 
(Molecular Dimensions) were used. Crystals of PL pi ° 
(ClllS)-Ub were observed under condition No. 9 of 
MD1-01, whereas no crystal of PL p o (WT)-Ub was ob¬ 
tained. Optimized crystals of PL p o (ClllS)-Ub were 
subsequently obtained within two days under the condi¬ 
tion: 22% w/v PEG 4000, 15% v/v 2-propanol, 0.1 mol/L 
tri-sodium citrate pH 4.8, and 10% glycerol. 2 pL of pro¬ 
tein and 2.5 pL of reservoir were mixed to equilibrate 
against 500 pL reservoir. 

Crystals were placed in a nitrogen-gas stream (100 K). 
A 3.16-A dataset was collected at wavelength 0.91841 A 
at beamline 14.2 of BESSY, Berlin (Mueller et al., 
2012). The program XDS (Kabsch, 2010) was used to 
process the diffraction data. The space group was found 
to be P6 3 , with unit-cell parameters a = b =138.14 A, c = 
57.59 A, y = 120°. Diffraction data statistics are shown 
in Table 1. 

Phase determination, model building and 
refinement 

The structure of the MERS-CoV PL pi ° (Cl 1 lS)-Ub com¬ 
plex was solved by molecular replacement using 
MOLREP (Vagin & Teplyakov, 2010). The program se¬ 
lected the MERS-CoV PL pro (Protein Data Bank (PDB) 
entry 4P16, Lei et al., 2014) as the first search model. 
Human ubiquitin (PDB entry: 1UBQ, Vijay-Kumar et 
al., 1987) was used as the second search model. The 
model of the complex was inspected and rebuilt using 
Coot (Emsley et al., 2010), and refined using phenix.re- 
fine (Headd et al., 2012). The refinement statistics are 


Table 1. Data collection and refinement statistics 


MERS-CoV PL pro (Cl 11 S)-Ub 


Data collection statistics 

Space group 

P6 3 

Unit-cell dimensions (A, °) 

a=£>=138.14, c=57.59 

y=120 

Wavelength (A) 

0.91841 

V m (A 3 /Da) 

3.59 

Solvent content (%) 

65.8 

Resolution range (A) 

45.22-3.16 (3.33-3.16) 

Number of unique reflections 

10955 (1577) 

R 1 
, 'pim 

0.050 (0.413) 

CC(1/2) 

0.997 (0.708) 

Completeness (%) 

99.9 (100) 

Mean I/a (1) 

12.8 (2.0) 

Multiplicity 

7.3 (7.5) 

Refinement statistics 

Rcryst (%) 2 

21.1 

Kfree (%f 

25.5 

No. of atoms 

Protein 

3054 

Ion 

1 

Ligand 

24 

Clashscore 2 

17 

r.m.s.deviation in bond lengths (A) 

0.012 

r.m.s.deviation in bond angles (°) 

1.945 

Average B-factor for all atoms (A 2 ) 

108 

Average B-factor for PL pro (A 2 ) 

95 

Average B-factor for Ub (A 2 ) 

145 

Ramachandran plot 

Residues in favored regions (%) 

95.30 

Residues in allowed regions (%) 

4.44 

Residues in outlier regions (%) 

0.26 


Notes: 1 Rpim (Weiss & Hilgenfeld, 1997). 2 R cryst = £ hk , 

| F 0 (hkl)-F c (hkl) | /£ hk) F 0 (hkl). R free was calculated for a test 
set of reflections (4.9%) omitted from the refinement. 

3 Cl-ash-score is defined as the number of clashes calculated 
for the model per 1000 atoms (including hydrogens) of the 
model (Chen et al., 2010). 

shown in Table 1 . The final model coordinates and struc¬ 
ture factors have been deposited in the PDB database 
with the code 4WUR. Figures (except for the supple¬ 
mental figure) have been prepared using Pymol (Schro- 
dinger; http://www.pymol.org/). 
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Kinetic assays of purified PL pro s 

All enzymatic assays were performed using a 96-well 
microtiter plate and the reaction buffer 20 mmol/L Tris- 
HC1, 150 mmol/L NaCl, pH 7.9, 2 mmol/L dithiothreitol 
(DTT). The fluorogenic substrates Cbz-Arg-Leu-Arg- 
Gly-Gly-7-amino-4-methylcoumarin (Z-RLRGG-AMC) 
(Bachem), Z-LRGG-AMC (BostonBiochem), and ubi- 
quitin-AMC (Ub-AMC) (BostonBiochem) were used. 
The fluorescence of free AMC with different concentra¬ 
tions (5 nmol/L-2.5 pmol/L) in reaction buffer was 
measured to generate a calibration curve, in order to con¬ 
vert the change of fluorescence intensity per unit of time, 
A(AFU) /s, into the amount of hydrolyzed substrate in 
pmol/L/s. 

The enzymatic cleavage reactions were run with an 
Flx800 fluorescence spectrophotometer (BioTek), to 
measure the increased fluorescence signal (k QX : 360 nm; 
X em : 460 nm) resulting from AMC release. Reactions 
were initiated by addition of the proteases to the reaction 
system. The peptide-hydrolysis kinetic assays were per¬ 
formed with the following conditions: 1 pmol/L MERS- 
CoV PL pro variant (D164E, D164A, D165E, F269Y), or 
10 pmol/L D165A, or 0.1 pmol/L SARS-CoV PL pro , with 
different concentrations (10, 20, 40, 80, 100 pmol/L) of 
Z-RLRGG-AMC or Z-LRGG-AMC in a final volume of 
100 pL at 25 °C. The kinetic curves for the proteases and 
their variants with the substrates Z-RLRGG-AMC or Z- 
LRGG-AMC were linear and the initial velocities also 
increased linearly with substrate concentration. No satur¬ 
ation was observed. Therefore, the data were fitted to the 
equation v/[E] tot = k app [ S], where k app approximates k ca /K M , 
as described previously (Barretto et al., 2005; Wojdyla et 
al., 2010). 

The deubiquitinating kinetic assays were performed 
under the following conditions: 0.1 pmol/L MERS-CoV 
PL pi ° wild-type or its variants D164E, D164A, D165E, 
F269Y, or 0.5 pmol/L D165A, or 0.025 pmol/L SARS- 
CoV PL pi ° were incubated with increasing concentra¬ 
tions (1, 2, 4, 8 pmol/L) of Ub-AMC in a final volume of 
50 pL, at 25 °C. Although PL pi ° actually cleaves the iso¬ 
peptide bond between the carboxyl group of the C-ter- 
minal Gly in Ub and the s-amino group of Lys in ubi- 
quitinated targets in vivo , we used the hydrolysis of Ub- 
AMC here to test the deubiquitinating activity in vitro. 
The kinetic curves of proteases and variants with the 
substrate Ub-AMC were hyperbolic and the initial velo¬ 
cities were not linear over the concentration of substrate. 
However, saturation was still not observed within a reas¬ 
onable time. Only when the ratio of protein to substrate 
was 1:1 or larger, were we able to achieve saturation 
within a limited time (data not shown). As the initial ve¬ 
locities did not increase in a strictly linear fashion with 
substrate concentration, application of the equation 
v/[E] tot = k [ S] to mimic k c JK M would lead to large 


standard errors. We were however able to fit the data to 
the Michaelis-Menten equation using the GraphPad 
Prism program (GraphPad Software), even though satur¬ 
ation could not be reached ( Supplemental Figure ). All 
assays were performed in duplicates. 

RESULTS 

Overall structure of MERS-CoV PL pro in complex 
with human ubiquitin 

The substrate-binding site of MERS-CoV PL pi ° features 
significant differences from those of the corresponding 
SARS-CoV enzyme and human ubiquitin-specific pro¬ 
teases (USPs, such as, USP14) (Hu et al., 2005; Chou et 
al., 2014; Ratia et al., 2014). It is therefore of interest to 
determine the crystal structure of the complex between 
MERS-CoV PL pi ° and its substrate, human ubiquitin. 
Hence, we crystallized the ubiquitin (Ub) complex of a 
MERS-CoV PL pi ° variant that had the active-site Cysl 11 
replaced by serine (Cl 1 IS) and determined the structure 
at 3.16 A (Figure 1A). 

There is one PL pi ° (Cl 1 lS)-Ub complex per asymmet¬ 
ric unit. Using the PDBePISA server (Krissinel & Hen- 
rick, 2007), the total interface region of MERS-CoV 
PL pi ° (Cl 1 lS)-Ub was determined as 813 A", close to the 
915 A 2 interface of SARS-CoV PL pi0 -Ubal (ubiquitin al¬ 
dehyde; PDB entry: 4MM3, Ratia et al., 2014) and the 
999 A 2 of SARS-CoV PL pro (Cl 12S)-Ub (PDB: 4M0W, 
Chou et al., 2014), but less than the 1503 A" observed for 
USP14-Ubal (PDB: 2AYO, Hu et al., 2005). The overall 
structure of MERS-CoV PL pi ° (Cl 1 IS) is very similar to 
that of the substrate-free PL pi ° (Lei et al., 2014), with a 
root-mean-square deviation (RMSD) of 0.91 A for cor¬ 
responding Ca atoms between these two structures. Nev¬ 
ertheless, two differences are immediately visible: 1), the 
zinc-finger motif has moved and closed in onto the Ub, 
compared to the free PL pi °. The zinc ion position has 
shifted by about 4 A and the largest deviation between 
the two structures is -6 A for Cys228 of the zinc-finger 
region (Figure IB); 2) the mobile loop 271 GIETAVG 77 , 
also named “BL2 loop”, is defined by clear main-chain 
electron density ( figure 1 ). This loop is disordered in 
substrate-free PL pi ° (Lei et al., 2014; Bailey-Elkin et al., 
2014). Compared to MERS-CoV PL pro , the BL2 loop 
267 GNYQCG" 72 is shorter by one residue in SARS-CoV 
PL pro . 

At the time when we deposited in the PDB the co¬ 
ordinates for the crystal structure of the non-covalent 
complex MERS-CoV (ClllS)-Ub (PDB: 4WUR), 
Bailey-Elkin et al. (2014) described two crystal struc¬ 
tures for a covalent complex MERS-CoV PL pi0 -Ub, in 
which an alkyl bromide group introduced at the C-ter- 
minus of Ub had formed a thioether with the active-site 
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Figure 1. Structure of MERS-CoV PL pro in complex with 
human ubiquitin (Ub). (A) Cartoon view of the overall 
complex structure. PL pro is shown in green and Ub in 
red. The N and C termini are marked by underlined let¬ 
ters, and the PL pro subdomains are divided by black 
dashed lines. The ubiquitin-like (Ubl) domain, thumb, 
palm, and fingers subdomains are indicated. The cata¬ 
lytic triad Cysl 11-His278-Asp293 is shown by spheres 
(yellow, blue, and red; Cys was replaced by Ser in this 
study). The zinc atom is displayed as a gray sphere. 
Five PL pro regions interacting with Ub are colored dark 
blue and labeled with Roman numbers (l-V). (B) Super¬ 
position of the fingers domain of Ub-bound MERS-CoV 
PL pro and the substrate-free enzyme (PDB: 4P16; Lei et 
al., 2014). The zinc atom is shifted by ~4 A, and the 
largest difference is ~6 A. (C) The BL2 loop (blue; 271 GI- 
ETAVG 277 ) of PL pro . A 2F 0 -F C electron density (gray; 1.0 
o) is displayed. The side-chains of residues in this loop 
are only partly defined by electron density. 


Cysl 11 of MERS-CoV PL pi °. These structures were in 
space groups P6 5 22 and P6 3 and were named “closed” 
and “open” PL pi0 -Ub complexes (PDB: 4RF0 and 4RF1), 
respectively (Bailey-Elkin et al., 2014). Even though the 
zinc-finger motif of PL pi ° in the former complex (space 
group P6 5 22) is closer to the Ub than in the P6 3 complex, 
both the “closed” and “open” complex show almost the 
same interactions between PL pi ° and Ub (Bailey-Elkin et 
al., 2014). The differences may be caused by crystal 
packing. Our non-covalent complex is similar to the 
“open” form of the covalent PL pi0 -Ub complex (RMSD = 
0.58 A for the PL pi ° and 0.85 A for Ub, based on all Ca 
atoms). 

All parts of the PL pi ° except for the ubiquitin-like 
(Ubl) domain interact with Ub. Most interactions in- 



Figure 2. Interactions between MERS-CoV PL pro and 
Ub. (A) Five regions of PL pro bind to two areas of Ub. 
The surface of PL pro is shown in wheat color. The five 
interacting regions are shown in brown, cyan, yellow, 
green, and purple, and they are also labeled by Roman 
numbers. The two interacting areas of Ub are marked 
“A” and “B”. Region A is highlighted by dark blue dots, 
and region B (the RLRGG sequence; Ub residues are in 
italics) is depicted in the ball-and-stick style. (B) Region 
A of Ub interacting with PL pro . A cartoon view of region 
A is shown in purple, and (33 and (34 are labeled. 
Residues of Ub (purple) and PL pro (green) are dis¬ 
played in the ball-and-stick style, and labeled in red and 
black, respectively. The 2F 0 -F C electron density (blue; 

1.0 a) of the side-chains of Arg42 and Aspl 65 is dis¬ 
played. Hydrogen bonds are indicated by black dashed 
lines, and the salt-bridge between R42 and D165 is de¬ 
picted by two red dashed lines. (C) The RLRGG bind¬ 
ing site of PL pro . The P5-P1 residues are shown in the 
ball-and-stick style. Regions I, II, IV, and V that interact 
with the P5-P1 residues are labeled. 

volve five surface regions of the PL pi ° and two regions of 
Ub ( figure 2A). These five regions of the PL pi ° are la¬ 
belled by Roman numbers: I, Leul06-Tyrl 12; II, 
Alai 62-Arg 16 8; III, Cys208-Val2 1 0; IV, 
Gly248-Pro250; V, Phe269-Tyr279. I and II are situ¬ 
ated in the thumb domain; III is in the fingers domain; 
and IV and V are in the palm domain ( figures 1A and 
2A). In the following and in the figure labels, residues of 
Ub are indicated in italics to distinguish them from 
residues of PL pi °. The two interacting regions of Ub are: 
A, Arg42-Gln49\ B, Arg72-Gly76 ( figure 2A). Region 
A of Ub consists of P3, P4, and the loop between p3 and 
P4 (Figure 2B). Regions II and III of PL pi ° interact with 
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region A of Ub. Region B comprises the five C-terminal 
residues, RLRGG , and is in contact with regions I, II, IV, 
and V of PL pi ° (Figure 2C). The C-terminal RLRGG mo¬ 
tif contributes the majority of the interactions with the 
PL pi °; the buried surface between Ub region B and PL pi ° 
is 477 A" (out of a total of 813 A 2 ). In order to make 
these interactions clear, we describe here in some detail 
the contacts between the PL pi ° and regions A and B of Ub. 

Interactions of MERS-CoV PL pro (ClIIS) with Ub 
region A 

Region A ( Arg42-Gln49 ) of Ub inserts into the space 
between the thumb and fingers domains of PL pi ° ( fig- 
ures 2A-B). Residues Ile44, Ala46, and Gly47 engage in 
hydrophobic interactions with Tyr209 and Val210 of re¬ 
gion III of PL pi ° ( figure 2B). The hydrophobic patch 
( [Ile44 , Ala46, and Gly47) is a common interaction re¬ 
gion utilized by Ub-binding proteins (Dikic et ah, 2009). 
In particular, the interaction of Ile44 with Val210 of 
PL pi ° is important for the deubiquitinase (DUB) but not 
for the protease activity. The variant V210R shows dra¬ 
matically reduced DUB activity, as demonstrated by 
Bailey-Elkin et al. (2014) (according to the numbering 
scheme of these authors, V210 is V1691). 

A salt-bridge exists between region A of Ub and the 
PL pi ° ( figure 2B), namely between the side-chains of 
Arg42 and Asp 165 (region II in PL pi °). In addition, a hy¬ 
drogen bond is formed between the main-chain O atom 
of Gly47 and the main-chain amide of Val210 (region 
III). In the SARS-CoV PL pro (Cl 12S)-Ub complex (Chou 
et al., 2014), Arg42 forms a salt-bridge with the negat¬ 
ively charged Glul68, a residue which is replaced by the 
positively charged Argl68 in MERS-CoV PL pi °. Con¬ 
sequently, the same salt-bridge cannot be formed in the 


MERS-CoV PL pro (C111S) -Ub complex. Instead, 
MERS-CoV PL pi ° has Asp 165 interacting with Arg42. In 
our deubiquitinating (DUB) kinetic assays, the D165A 
variant shows a dramatically reduced DUB activity; the 
k cat IK M * s a bout 78-fold decreased compared to that of 
the wild-type (Table ). The Km value of the D165A 
variant is about 4-fold higher than that of the wild-type 
MERS-CoV PL pi °, suggesting that this amino-acid re¬ 
placement reduces the Ub binding affinity. Meanwhile, 
the D165E amino-acid replacement shows a catalytic ef¬ 
ficiency comparable to wild-type towards Ub-7-amino-4- 
methylcoumarin (Ub-AMC) ( Table 1 ). However, we no¬ 
ticed that the K M value of D165E for the DUB activity is 
about 2-fold larger than for the WT enzyme. Although 
Glul65 can mimic Asp 165 here, we propose that the 
longer side-chain of Glu may fit less perfectly compared 
to Asp. These results indicate that the salt-bridge 
between Asp 165 and Arg42 could be important for the 
PL pi0 ’s DUB activity, in addition to the interaction 
between Asp 165 and the P4-amide group (see below). 

Interactions of MERS-CoV PL pro (ClIIS) with Ub 
region B 

Region B of Ub comprises the five C-terminal residues, 
RLRGG. These five residues bind to the narrow active- 
site channel between the thumb and palm domains of the 
PL pi ° ( figure 2C ). They mainly interact with regions I, II, 
IV, and V of the protease. Residues RLRGG are compat¬ 
ible with the PL pi ° cleavage motif, (R/K) (L/I) XGG 
(P5-P1), in the MERS-CoV polyproteins; their interac¬ 
tions with the PL pi ° are discussed here in terms of sub¬ 
sites SI to S5. 

SI and S2 subsites. Pro 163 (region II of PL pi °) and the 


Table 2. Kinetic parameters of MERS-CoV PL pro and SARS-CoV PL pro 



Z-RLRGG-AMC 

Z-LRGG-AMC 


Ub-AMC 

k app (Mmol -1 min -1 )* 

K m ((jmol/L) 

kcat ( min 1 ) 

k J K M (Mmol/L -1 min -1 ) 

MERS-CoV PL pro 

WT*' 1 

(1.3 ± 0.1) x 10 -3 

(1.0 ± 0.01) x 10 -3 

6.1 ±0.7 

8.4 ±0.7 

1.4 ±0.2 

F269Y 

(2.0 ± 0.1) x 10 -3 

(1.4 ± 0.2) x 10 -3 

8.5 ±0.1 

22.4 ± 1.5 

2.6 ±0.2 

D165E 

(1.1 ± 0.2) x 10 -3 

(8.4 ± 1.1) x 10 -4 

15.9 ±2.7 

27.7 ± 1.2 

1.7 ±0.3 

D165A 

(8.2 ± 1.6) x 10 -5 

(1.0 ± 0.2) x 10 -4 

23.9 ±6.2 

0.42 ± 0.04 

0.018 ±0.005 

D164E 

(1.5 ± 0.1) x 10 -3 

(8.6 ± 1.8) x 10 -4 

15.7 ±4.4 

15.3 ±3.3 

1.0 ±0.3 

D164A 

(2.9 ± 0.1) x 10 -4 

(6.5 ± 0.8) x 10 -4 

6.5 ±0.6 

2.5 ±0.01 

0.40 ± 0.04 

WT # ’ 2 

(3.0 ± 0.1) x 10 -3 


14.3 ±2.0 

18.8 ± 1.2 

1.3 ±0.2 

SARS-CoV PL pro 

WT*' 1 

0.23 ± 0.04 

(3.5 ± 0.5) x 10 -2 

57.2 ± 10.9 

197.4 ±28.0 

3.5 ±0.8 

WT # ’ 2 

0.3 ±0.1 


50.6 ± 7.4 

75.9 ±8.1 

1.5 ±0.3 

Notes: # : Wild type; 

*: k app is an approximation to k ca /K M ; b 

Our results (this work); 2 : Baez-Santos et al., 2014b. 
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side-chains of Asnl09 and Tyrl 12 (located in region I) 
form a space-restricted SI site to accommodate Gly76 
(PI). The carbonyl oxygen atom of Gly277 (region V) 
accepts a hydrogen bond from the amide of Gly76 ( fig¬ 
ure 3). 

The side-chains of Tyrl 12 (region I) and Phe269, 
Val276 as well as Tyr279 (region V), and Gly277 form a 
restricted space for Gly75 (P2). Two hydrogen bonds 
link the main-chain at Asp 164 of the PL pi ° and Gly75 
( ugure 3). The SI and S2 sites of the protease are too 
small to accommodate any other residue but glycine. 

S3 subsite. The main-chain O atom of Arg74 (P3) ac¬ 
cepts a hydrogen bond from the Gly277 amide. In the 
complexes SARS-CoV PL pro (Cl 12S)-Ub or SARS-CoV 
PL pi0 -Ubal, the main chain at Arg74 forms two hydrogen 
bonds (Chou et al., 2014; Ratia et al., 2014), namely with 
the amide of Gly272 and with the hydroxyl group of 
Tyr265. The former hydrogen bond is conserved in the 
MERS-CoV PL pi ° (Cl 1 lS)-Ub complex, but the latter is 
not. Tyr265 of SARS-CoV is replaced by Phe269 in 
MERS-CoV, which lacks the ability to form a hydrogen 
bond with the main-chain amide of Arg74. In agreement 
with this difference, the DUB activity of SARS-CoV 



Figure 3. Details of the interactions between the C-ter- 
minal RLRGG of Ub and MERS-CoV PL pro . RLRGG 
residues are shown in purple in the ball-and-stick style, 
and they are labeled in red. For clarity, the 2F 0 -F C elec¬ 
tron density (blue; 1.0 a) of the RLRGG main chain is 
shown. PL pro residues are displayed in green in the ball- 
and-stick style, and they are labeled in black. Hydrogen 
bonds are displayed as black dashed lines, and the 
salt-bridge between R72 and D164 is depicted as two 
red dashed lines. 


PL pi ° is 2.5-fold higher than that of MERS-CoV PL pi ° in 
our kinetic assay; furthermore, the MERS-CoV PL pi ° 
F269Y amino-acid replacement leads to enhancements 
by about 1.5-, 1.4-, and 1.9-fold of the hydrolytic activit¬ 
ies towards carbobenzoxy-Arg-Leu-Arg-Gly-Gly-7- 
amino-4-methylcoumarin (Z-RLRGG-AMC), Z-LRGG- 
AMC, and Ub-AMC, respectively ( Table T ). 

The side-chain of Arg74 is exposed to the solvent in 
our MERS-CoV PL pi ° (Cl 1 lS)-Ub complex. This side- 
chain shows remarkable variability in the interactions it 
makes in the different complexes. In the “open” but not 
in the “closed” covalent MERS-CoV PL pi0 -Ub complex 
(Bailey-Elkin et al., 2014), it donates a hydrogen bond to 
the main-chain carbonyl oxygen of Thrl755 of the BL2 
loop (corresponding to Thr274 in our numbering 
scheme). In SARS-CoV PL pi0 -Ubal (Ratia et al., 2014) 
but not in the SARS-CoV PL pi ° (Cl 12S)-Ub complex 
(Chou et al., 2014), the side-chain of Arg74 forms a hy¬ 
drogen bond with the main-chain carbonyl oxygen of 
Gln270. Instead, Arg74 is involved in a relatively weak 
salt-bridge with Glul62 in the SARS-CoV PL pi ° 
(Cl 12S)-Ub complex (Chou et al., 2014). None of these 
interactions exist in our MERS-CoV PL pi0 -Ub complex. 

54 subsite. The main-chain amide of Leu73 (P4) donates 
a hydrogen bond to the side-chain of Asp 165 (region II) 
( ugure 3). This hydrogen bond is conserved in all the 
complexes compared here. In our peptide-cleavage as¬ 
say, the D165A variant shows about 16-fold and 10-fold 
lower activities towards substrates Z-RLRGG-AMC and 
Z-LRGG-AMC, respectively, indicating that Asp 165 is 
not only important for interacting with Arg42 of Ub (see 
above). 

The side-chain of Leu73 is embedded in a hydrophobic 
pocket which is formed by the Cp atom of Asp 165, the 
side-chain of Pro250 (region IV), Phe269, as well as the 
Cp and Cy atoms of Glu273 (region V). Asp 165 and 
Pro250 are conserved in SARS-CoV PL pi ° (Asp 165 and 
Pro249). Phe269 is replaced by Tyr265, and Glu273, 
situated in the BL2 loop, is replaced by Tyr269. 
However, the side-chains of all these non-conserved 
residues possess the ability to provide a hydrophobic en¬ 
vironment to accommodate Leu73. 

55 subsite. The side-chain of Arg72 (P5) is located 
between the PL pi ° thumb domain and region A of ubi- 
quitin. It forms a salt-bridge with the side-chain of 
Asp 164 in MERS-CoV PL pi °. In addition, the guanidin- 
um group of Arg72 may be involved in a 7T-7T interaction 
with that of Argl68 ( ? igure ). These interactions have 
also been described for the covalent complex of MERS- 
CoV PL pi ° with Ub (Bailey-Elkin et al., 2014). Arg72 is 
not subject to strict space limitations; accordingly, this 
residue displays different binding patterns in the two 
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SARS-CoV PL pi0 -Ub complexes. In the covalent SARS- 
CoV PL pi0 -Ubal complex (Ratia et al., 2014), Arg72 
forms a salt-bridge with Glul68. In the non-covalent 
SARS-CoV PL pro (Cl 12S)-Ub complex (Chou et al., 
2014), Arg72 is exposed to the solvent and does not in¬ 
teract with Glul68 (Arg42 instead forms a salt-bridge 
with Glul68, as mentioned above). In our kinetic assay, 
the D164A variant of MERS-CoV PL pi ° displays a -4.5- 
fold and a -3.5-fold reduced activity, respectively, for Z- 
RLRGG-AMC and Ub-AMC (Table 2). For Z-LRGG- 
AMC, the activity is decreased just a little (by about 
-1.5-fold), because there is no P5-Arg in this substrate 
( Table T ). These data demonstrate that Asp 164 is import¬ 
ant for the interaction with Arg72. 

In summary, the main-chain heteroatoms of P5-P1 
form a hydrogen-bonding network with PL pi °. The bind¬ 
ing characteristics of PI, P2, and P4 are conserved in all 
MERS-CoV and SARS-CoV PL pro -Ub complexes. 
However, P3 -Arg and P5 -Arg assume binding patterns 
that differ between the various MERS-CoV and SARS- 
CoV PL pro -Ub complexes. 

DISCUSSION 

Viral proteins are likely to possess multiple functions, as 
exemplified by non-structural protein 1 (NS1) of influ¬ 
enza A viruses (Hale et al., 2008), the nucleocapsid (N) 
protein of coronaviruses (Chang et al., 2014), or the 
Nspl4 exonuclease-guanyl-7-methyltransferase of 
coronaviruses (Chen et al., 2009). Exhibiting DUB and 
proteolytic activities, the papain-like proteases of 
coronaviruses are no exception here. Although the over¬ 
all folds are conserved, the enzyme activity and sub¬ 
strate-binding modes of CoV PL pio s differ in detail. 
Therefore, no coronavirus PL pi ° can be considered a gen¬ 
eral model for all its homologues (Baez-Santos et al., 
2014b). 

MERS-CoV PL pro is not optimized for catalytic 
efficiency 

We have previously noticed that the oxyanion hole of the 
MERS-CoV PL pi ° appears to be deficient (Lei et al., 
2014; also see the discussion below). Similarly, the re¬ 
cognition of ubiquitin by the enzyme appears to be sub- 
optimal. Thus, the main-chain amide of the P3 -Arg 
residue has no hydrogen-bonding partner on the MERS- 
CoV PL pro , because the near-by side-chain of Phe269 is 
incapable of accepting an H-bond. The corresponding 
residue is Tyr265 in the SARS-CoV PL pro , which is per¬ 
fectly positioned to accept the hydrogen bond from the 
P3-amide. Indeed, when we replaced Phe269 by Tyr in 
MERS-CoV PL pro , the DUB activity of the enzyme in¬ 
creased by a factor of almost 2 and the peptidolytic activ¬ 


ity by -1.5. The evolution of viral enzymes is not neces¬ 
sarily driven by optimization of catalytic efficiency. This 
is particularly true for viral proteases that have to ensure 
the availability of non-structural proteins in the correct 
temporal order when they cleave them out of the viral 
polyproteins; in fact, too rapid a polyprotein processing 
might be counterproductive. On the other hand, a more 
efficient DUB activity should help the virus in counter¬ 
acting the innate immune response of the host cell. As 
the binding of the P5-P1 residues of ubiquitin and of the 
polyprotein cleavage sites obviously influences both the 
DUB and proteolytic activities of the PL pi °, the subop- 
timal catalytic efficiency that we observe may be a con¬ 
sequence of a compromise between the requirements of 
the two activities. 

With regard to the oxyanion hole, we previously pro¬ 
posed that the backbone amide of Asnl09 (located in a 
P-turn connecting P7 and a4, 7 igure 4) may contribute to 
the stabilization of the oxyanion intermediate in PL pi ° 
catalysis, along with the main-chain amide of the active 
site Cysl 11, although this may require a slight rearrange¬ 
ment of this P-turn (Lei et al., 2014). In our non-cova¬ 
lent MERS-CoV PL pi0 -Ub complex, we do not see any 
rearrangement of the P-turn. More or less in agreement 



Figure 4. The side-chain of Asn109 is unlikely to con¬ 
tribute to the oxyanion hole of MERS-CoV PL pro . PL pro is 
shown in green in cartoon view. N and C termini are 
marked by underlined letters. Strand P7, helix a4, and 
the loop between them are colored in blue. Asn109 and 
Gly 161 are indicated in ball-and-stick style in cyan and 
purple, respectively, and they are labeled in black. The 
hydrogen bond between them is displayed as a black 
dashed line. The active-site Cysl 11 (replaced by Ser in 
this study) is shown in green in ball-and-stick style. 
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with our suggestion, Bailey-Elkin et al. (2014) proposed 
on the basis of their covalent MERS-CoV PL pi0 -Ub com¬ 
plex structure that the main-chain amides of Asnl590, 
Asnl591, and Cysl592 (corresponding to Asnl09, 
Asnl 10, and Cysl 11 in our numbering scheme) form the 
oxyanion hole. On the other hand, Lee et al. (2015) ar¬ 
gued that the side-chain of Asnl 09 could contribute to 
the oxyanion hole. They found that the N109A replace¬ 
ment completely abolished enzyme activity. As we re¬ 
ported earlier (Lei et al., 2014), the side-chain amide of 
Asnl09 makes a strong hydrogen bond with the con¬ 
served Glyl61 ([figure 4). Any reorientation of the 
Asnl09 side-chain towards the oxyanion would require a 
disruption of this strong interaction; this is not very 
likely. In conclusion, lacking a side-chain in the proper 
spatial orientation and capable of donating a hydrogen 
bond to the oxyanion (such as Trpl07 in SARS-CoV 
PL pi °), the oxyanion hole of MERS-CoV PL pi ° seems to 
be deficient. 

Unique features of Ub recognition by MERS- 
CoV PL pro 

Apart from the less than optimum binding of the P3 
residue, there are other differences in the way MERS- 
CoV PL pi ° and SARS-CoV PL pi ° recognize human ubi- 
quitin. The formation of a salt-bridge between Arg72 , the 
P5 side-chain of Ub, and Asp 164 is unique for MERS- 
CoV PL pi °. The same interaction exists in the two cova¬ 
lent MERS-CoV PL pi0 -Ub complexes (Bailey-Elkin et 
al., 2014). This binding mode is very different from the 
Glut68 - Arg72 salt-bridge in the SARS-CoV PL pro -Ubal 
complex (Ratia et al., 2014). Glul68 is conserved in 
HCoV-NL63 PL2 pi ° and replaced by Asp in the HCoV- 
229E, HCoV-OC43, and HCoV-HKUl PL2 pr °s (for se¬ 
quence alignments, see Barretto et al., 2005; Baez-San- 
tos et al., 2014b), so the same type of interaction is likely 
to be realized in the Ub complexes of these enzymes. 
However, the corresponding residue in MERS-CoV PL pi ° 
is Argl68; hence, Asp 164 is used instead to bind Arg72. 
Asp 164 is in fact unique in MERS-CoV. It is replaced by 
Gly in SARS-CoV PL pro and the PL2 pr °s of HCoV NL63 
and HCoV 229E, by Ala in HCoV-OC43 PL2 pro , and by 
Ser in HCoV-HKUl PL2 pi °. Our kinetic results for the 
D164A replacement (see Results) emphasize the import¬ 
ance of the unique Asp 164 residue. 

The role(s) of the ubiquitin-like (Ubl) domain 

In our PL pi ° (Cl 1 lS)-Ub complex, the Ubl domain of 
PL pi ° shows no interaction with ubiquitin. The relative 
orientation of the Ubl domain in the substrate-bound 
PL pi ° is the same as in substrate-free PL pi °. The Ubl do¬ 
main of MERS-CoV PL pi ° is not required for the IFN 
antagonism activities (Baez-Santos et al., 2014b), but it 
is required in case of SARS-CoV PL pi ° (Frieman et al., 


2009). Recently, Mielech et al. (2015) reported that the 
Ubl domain of mouse hepatitis virus (MHV) is an im¬ 
portant modulator of PL pi ° stability and viral pathogenes¬ 
is. Although the Ubl domain shows variable effects in 
different Co Vs, the high degree of conservation of the 
domain throughout the family suggests that it may play a 
common biological role. One possible function is that the 
Ubl might be involved in protein-protein interactions. 
Ubiquitin-like domains are known to function as binding 
modules in such interactions. For example, the kinase 
Raf contains a Ubl domain for interaction with human 
Ras (Fetics et al., 2015). Also, human ubiquitin-specific 
protease (USP) 7 includes five Ubl (1-5) domains, of 
which the second is bound by the Herpes simplex virus-1 
immediate-early protein ICP0 to counteract the intrinsic 
antiviral response of the host cell (Pfoh et al., 2015). 

PL pro inhibitors 

The BL2 loop of SARS-CoV PL pi ° is important for bind¬ 
ing inhibitors (Baez-Santos et al., 2014a; Lee et al., 
2015). This loop is variable in different CoV PL pi °s. A 
potent inhibitor of the SARS-CoV PL pro , N-[(3-fluoro- 
phenyl) methyl]-1 -[(1 R) -1-naphthalen-1-ylethyl] piperid- 
ine-4-carboxamide (compound 3k, IC 50 = 0.15 ± 0.01 
pmol/L) was found to have no effect on the MERS-CoV 
enzyme (Baez-Santos et al., 2014b). Tyr269 and Gln270 
of SARS-CoV PL pi °, which are important for binding this 
inhibitor (Baez-Santos et al., 2014a), are replaced by 
Glu273 and Ala275 in the MERS-CoV protease. It seems 
that this structural difference in the BL2 loop has a re¬ 
markable impact on the effectiveness of the inhibitor. 
The structure of the MERS-CoV PL pi ° (Cl 1 lS)-Ub com¬ 
plex presented here will facilitate virtual screening of 
chemical libraries for specific anti - MERS-CoV PL pi ° 
inhibitors (Hilgenfeld, 2014). 

The di-Ub site of PL pro 

SARS-CoV and MERS-CoV PL pro s can digest K48- and 
K63-linked (poly)ubiquitin chains and remove ISG15 
from proteins covalently linked to it (Ratia et al., 2014, 
Baez-Santos et al., 2014b). SARS-CoV PL pi ° prefers 
binding of K48-Ub 2 and ISG15 over mono-ubiquitin 
(Ratia et al., 2014). All evidence suggests that on the 
PL pro , at least two ubiquitin-binding sites (or a binding 
site for diubiquitin-like molecules such as ISG15) exist. 
Ratia et al. (2014) proposed hypothetic models for com¬ 
plexes of SARS-CoV PL pro with K48-Ub 2 and ISG15. 
These authors identified two major hydrophobic binding 
sites on the PL pi ° for the first Ub (Ubl) and the second 
Ub (Ub2). The binding site for Ubl comprises Met209, 
Pro248, and Pro249. This hydrophobic patch is con¬ 
served in MERS-CoV PL pi °, although the residues are 
not exactly the same. Met209 of SARS-CoV PL pi ° is re¬ 
placed by Val210 (region III) in MERS-CoV. Pro248 is 


www.virosin.org 


9 



Structure of MERS-CoV PL pro -ubiquitin complex 


replaced by Thr249, whereas Pro249 is conserved 
(Pro250; region V) (figure 5 ). The hypothetic second 
binding site for Ub2 (also named “ridge” region, Ratia et 
ah, 2014) is located to the first a helix (a2) in the thumb 
domain, including residues Phe70, His74, and Leu76 of 
SARS-CoV PL pi °. According to a structural alignment of 
the SARS-CoV and MERS-CoV PL pr °s, Phe70, His74, 
and Leu76 are changed to Lys69, Gly73, and Val75, re¬ 
spectively ( figure f ). In SARS-CoV PL pro , the F70S and 
F70A replacements lost the affinity to K48-Ub 2 and 
ISG15 in vitro (Ratia et ah, 2014); therefore, the pres¬ 
ence of Fys69 in MERS-CoV PF pi ° instead of Phe70 
strongly suggests that this enzyme should bind Ub2 in a 
different way. Bailey-Elkin et al. (2014) predicted 
Asnl673 and Vall674 of the fingers subdomain (corres¬ 
ponding to Asnl92 and Vail93 in our PF pi °) as the distal 
Ub site of K63 di-Ub, but they found that their DUB 
activity data do not support their prediction for the Ub2 
binding site. The structure of the MERS-CoV PL pi0 -Ub 
complex reported here reveals the exact Ubl binding site 
on the PL pro , but the Ub2 binding site should be identi¬ 
fied by crystallizing the PL pi ° in complex with poly- or 
di-ubiquitin. 

In summary, the crystal structure of the MERS-CoV 
PL pi0 -Ub complex provides valuable information that 
helps understand the multiple functions of coronavirus 
papain-like proteases. Mutational studies additionally 
highlight features of the MERS-CoV PL pi °. The different 
substrate-binding patterns should be kept in mind when 
designing inhibitors for PL pi °s of Co Vs, even though the 
overall structures of these enzymes are conserved. Fur¬ 
thermore, it would be helpful to obtain the structure of 



Figure 5. The two Ub-binding sites of PL pro . The Ubl 
and the proposed Ub2 (according to a structural align¬ 
ment with the SARS-CoV PL pro ; Ratia et al., 2014) bind¬ 
ing sites are depicted as purple dots. The N and C ter¬ 
mini of PL pro are marked by underlined letters. All 
residues related to the two Ub binding sites are labeled. 
The catalytic triad Cysl 11-His278-Asp293 is indicated 
by yellow, blue, and red spheres. 


PL pi ° in complex with di-Ub or ISG15 in the future. 
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